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ABSTRACT 
> 

O . 

Q , Using a novel two-dimensional coordinate system, we have derived a 

^ ! particularly simple way to express the redshift distortions in galaxy redshift 

t^^ ! surveys with arbitrary geometry in closed form. This method provides an almost 

- ' ideal way to measure the value of /? = f^Q^/fo in wide area surveys, since all 



o 



pairs in the survey can be used for the analysis. In the limit of small angles, 
this result straightforwardly reduces to the plane-parallel approximation. This 
expansion can also be used together with more sophisticated methods such 

j^ \ as for the calculation of Karhunen-Loeve eigenvectors in redshift space for an 

arbitrary survey geometry. Therefore, these results should provide for more 
k> \ precise methods in which to measure the large scale power spectrum and the 

H \ value of /3 simultaneously. 



Subject headings: cosmology: large-scale structure of the universe 



1. Introduction 

The fact that the correlation function in redshift space is distorted from the real space 
correlations is well known. Davis and Peebles (1983) have demonstrated quantitatively 
that the correlation function ,^(7r,cr), expressed as a function of the line-of-sight (vr) and 
perpendicular [a) separations, has a significant anisotropy. By measuring the elongation of 
this function along the (vr) axis, a value of the one-dimensional velocity dispersion can be 
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inferred. They found that the velocity distribution was well described by an exponential 
distribution and that the value of the dispersion was relatively independent of scale, at 
around 350 km/s. Peebles (1980) gives a nice summary why in linear theory the typical 
velocities of the galaxies scale with fig'^. With the advent of biased galaxy formation 
(Bardeen et al 1986), the velocities have been rescaled by a bias factor 6, leading to velocities 
dependent upon the combination (3 = VlP'^ /h. 

In Kaiser's (1987) pioneering work, the effect of the infall due to linear theory was 
identified and it was shown that in the limit that the lines of sight to the two galaxies are 
approximately parallel to each other, the power spectrum P{k) is enhanced as a function of 
the directional cosine /i between the wave vector and the line-of-sight as 

P\k,ii) = P{k){l + Pl^'f- (1) 

In related work, Lilje and Efstathiou(1989) calculated the angular average of the 
redshift correlation function directly related to this expression. IIamilton(1992) expanded 
the redshift space correlation function into components, multiplied with the angular 
multipoles. He showed that only the quadrupole and hexadecapole terms arise in linear 
theory, in agreement with Kaiser's expression. Other promising approaches have involved 
an expansion into orthogonal eigenfunctions (Heavens and Taylor 1994), or restricting the 
surveys to small opening angles so that the plane-parallel approximation still holds (Cole, 
Fisher and Weinberg 1994,1995). These papers attempted to measure the parameter (3 from 
the quadrupole to monopole ratio. 

Zaroubi and Hoffmann (1996) outlined how to compute a linear expansion of the 
redshift space correlation function for a general geometrical configuration and provided 
numerical estimates of the redshift distortions. Additionally, Hamilton and Culhane (1996), 
hereafter HC96, have introduced a novel integral transform, rotationally invariant and 
commuting with the redshift distortion operator, which makes the transformation very 
elegant and simple. 

Here we present simple and intuitive expressions to calculate the redshift distortions 
which exploit the fact that the inherent geometry of the problem is two dimensional. 
Indeed, by rotating all pairs in a redshift survey to a common plane and analyzing the 
redshift distortions in that plane, we obtain particularly simple analytic results. These 
results smoothly approach the Kaiser/Hamilton plane-parallel limit in which angle between 
the galaxy pair in a survey is small. Additionally, these expressions allow for a much more 
productive use of redshift surveys since all galaxy pairs can be utilized in a given wide angle 
survey. It can also be used in adjunct with more sophisticated methods of analysis such as 
construction of eigenvectors for a KL analysis (Vogeley and Szalay 1996). 
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2. The Correlation Function in Redshift Space 

2.1. The Coordinate System 

The symmetries in the geometric configuration of redshift surveys have been elegantly 
described in HC96, who recognized that the correlation function in redshift space is tied to 
a triangle formed by the observer and the two galaxies which lie upon two arbitrary lines of 
sight. The correlation function should be a function of this triangle, but otherwise invariant 
under rotations about the observer. Such invariant triangles can be characterized by one 
size and two shape parameters. We will use two angles to describe the shape of the triangle, 
and express the scale dependence of the correlation function as a function of r = |r|, the 
distance between the two points. 

In this coordinate system, which is illustrated in Fig. 1, the two points are located at 
points ri = rifi and r2 = T2r2. The observer is at the origin, at rest in comoving frame. The 
angle between the two normal vectors is 29, that is rir2 = cos 2^. The separation between 
the two points is r = ri — r2, with a normal vector r. The angles between ri, r2 and r are 
7i and 72 respectively. The symmetry axis for ri and r2 is found by halving the angle 29. 
The angle between this symmetry axis and r is 7. We assume that ri > r2, thus 72 > 71. 
The angles are related to each other as 71 + 72 = 27, and 72 — 71 = 29. This introduction 
implies our choice of the two angles, 7 and 9., which together with |r| completely describe 
the shape of the triangle. These angles are particularly convenient since 9 is given by our 
geometry and 7 smoothly approaches the angle familiar in the plane-parallel limit as 9 
tends to 0. There is another geometric relation between the two angles: if 9 is not 0, then 
the smallest value 7 can take is 9, and its largest value is vr — ^. This is not due to the 
choice of coordinates but rather to the geometry of the problem. 

We will provide an expression for the distortions of the correlation function as a series 
expansion in terms of the angle 7, where the coefficients have a simple r and 9 dependence, 
together with a dependence on the parameter (3 that we wish to measure. So far our 
approach has been quite similar to most earlier work except for the introduction of a 
convenient set of angles which will simplify the results considerably. 

It has been a persistent problem to get a simple power spectrum or correlation function 
beyond the Kaiser (1987) plane-parallel limit. The difficulty lies in that the Fourier space 
transfer function, the ratio of the distorted and undistorted power spectrum contains a 
strong mode-coupling (Zaroubi and Hoffmann 1996). Thus, the transfer function is not 
multiplicative but contains an integral over a non-local kernel with a strong /c- dependence. 
This dependence strongly affects the multipole coefficients which would otherwise make 
a precise determination of (3 quite straightforward. One alternative has been to use only 
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galaxy pairs at small angular separations (e.g. Cole et al 1994, 1995). 

It is interesting to consider what the reasons are for this non-locality. With a 
finite 9, the original symmetry of the problem is lost, and the geometry is inherently 
two-dimensional. This effect is thus arising from 'aliasing' (Kaiser and Peacock 1991, Szalay 
et al 1993, Landy et al 1996), i.e. projections of the spherically isotropic power to the lower 
dimensions of the survey geometry. Interestingly enough, even though the power spectrum 
is definitely non-local, the redshift space correlation function can still be computed in 
spherical coordinates, and it can be expressed in a closed form, as will be shown below. 

In the following Section 2.2, we will first derive a representation of the redshift 
distortion problem in terms of spherical tensors. This development, being a representation 
independent expansion, is independent of our choice of coordinate system described above. 
In Section 2.3, we will reintroduce our coordinates and show how the coefficients of this 
expansion can be most economically expressed using these coordinates. 



2.2. Expansion in Spherical Tensors 

We can assume without a loss of generality, that the observer is at rest with respect 
to the CMB frame. As shown by HC96, this effect can be trivially included. The linear 
expansion of the overdensity at a redshift space coordinate s relates to the real-space 
overdensity at r, with s = r -|- -ur, as a function of the radial velocity u and its line-of-sight 
derivative du/dr, as first introduced by Kaiser (1987), 

^W(r)=<5M(r)_^!fM_tt(r)M(r), (2) 

or 

Here a{r) = {2 + (9 In 0(r)/(9 In r)/r, where 0(r) is the selection function which is a slowly 
varying function of r, the radial distance from the observer. Due to the fact that the 
velocity scale is much smaller than the typical depth of todays redshift surveys, this term 
is very small and is generally ignored. We will include this term for completeness in our 
results but would like to make one other point why this term is additionally small. If the 
redshift survey is defined by a boundary on the sky, a given selection function, and a range 
of the radial coordinate (a, 6), then a will be averaged, weighted by the selection function. 
Integrating by parts one can show that this term is small and becomes zero as the lower 
limit becomes and the upper limit oo, the case of full surveys. 

The peculiar velocity can be written concisely as a Fourier integral, with the usual 
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(3 = n^/h, 



as 



Mr 



i(3 



d^k kf 



— e*'^'"5(k). 



(27r)3 k 
In order to simplify subsequent calculations, we introduce here the spherical tensor A as 

(i^k 



(3) 



^r(r) 



{2nf 



(iA;)-"Pi(kr)e*'^'-5(k) 



One can then reexpress the redshift-space overdensity at the point r in terms of A as 
5W(r) = (iH- f ) Aliv) + ^-pAl{v)+Pa{r)A\{v). 



(4) 



(5) 



This use of spherical coordinates makes it easy to express the correlation of the different A 
components, at positions ri and Y2 as 

(fk 



(,^)-{n,+n.)p^^(kfi)P^,(kr2)e^'^^P(A;). (6) 



where P^ is the Legendre polynomial, r = ri — r2 and P{k) is the power spectrum with 
the definition (5(k)5(k')) = {27r)^5^(k + \i.')P{k). The redshift-space two-point correlation 
function of finite-angle is given in terms of SJ^i^, where we will use a, = a{ri). 



^^'\r^,T,] 



(1 + f ) ' ^o°o + l^'S',, + ^/5 (1 + f ) {s'o2 + SI 



+ ^ ( ^ + ^ ) '^I'^io + "S'S'oi + -/3^ ai5'J2 + «25'2i - /3^aia25'n (7) 

Both the Legendre polynomials and the plane wave e^^'^ can be expanded in terms of 
spherical harmonics as 



m=—l 
00 I 

e^''" = 47rE E ^'jiikr)Y,Um^{i), 

Z=0 m=-i 



(9) 



where ji and YJm, are spherical Bessel functions and spherical harmonics, respectively. We 
can express the integral of three l^m's over (ii7k with the Wigner 3j-symbols. We also 
introduce the bipolar spherical harmonics X;^^(fi,r2) which transform as a spherical 
harmonic with L, M with respect to rotations (Varshalovich et al 1988). 



^zt"(ri>r2) 



-1 



di-h-M 



V^Sl+T E (^ !! \. ]Y,,rrAri)Y,,mAh) (10) 

mi,m2 \ 
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With Ylm{^) these can form a rotationally invariant scalar, dependent on L, li, I2, where A 
indicates the shape of the triangle formed by the three unit vectors. Any quantity that is a 
scalar function of A is independent of what angles or parameters we may choose to describe 
the triangle shape. This invariant scalar B is symmetric in /i,/2- 



"^"'"' ^/,2,.L,.n (ooo)g^^?•'-^--)^^"'^)■ 



(11) 



^(2/i + l)(2/2 + l) 
With this notation, equation (|) reduces to 

^y. = (4vr)^/^(-l)'^ E^"""^y.(A)ef (0 (12) 



where 



^t\r) = ^Jdkek--jL{kr)P{k), (13) 



2.3. Expansion of the Correlation Function 

The above expression shows the correlation function can be written in terms of a 
series, factorized into size r and shape A dependent terms, providing a representation 
independent expansion. Equations (^ and (|1^) give an expression of the redshift-space 
correlation function. Re-expanding in terms of Cl (r), the A dependence is contained in 
the coefficients Cjj(A): 

e(^) (r. A) = cooeS°^ + Co2d°^ + Co,^P + ciid'^ + ciad'^ + C2od'^ + ^22^^ , (14) 

It is at this point that the use of any particular coordinate system becomes paramount. 
Since the coefficients Cjj(A) are functions of the shape of the triangle under consideration, 
a judicious use of coordinates can simplify the problem considerably, as will be shown 
below. The terms proportional to sin^(6') have been saparated, since these disappear in the 
plane-parallel limit. After a somewhat tedious reduction procedure, most of it carried out 
in Mathematica, we obtain the following results: 

Coo = 1 + ^/5 + ^P' - ^/5' cos^ 9 sin2 9, (15) 

C02 = -{p+ ^P') cos 20P2(cos7) -l{^-\^' + \^' sin^ 9^ sin^ 9, (16) 

Co4 = ^/3^P4(cos7) - ^(3^ sin^ 9P,{cos^) - ^-(3' ( A _ ^ ^.^2 ^^ ^^^2 q ^^7) 
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and 

cii = {ai-a2) (p + -P"^ - -P'^sm'^9jcos9Pi{cos'j) 

+ {ai + a2)(p + -p^--p^cos^9]sm9Pi{sm-f) (18) 

Ci3 = -(ai — a2)/3^cos^[sin^ ^Pi(cos7) — 2P3(cos7)] 



+ -(ai + a2)/5^sin^[cos2^Pi(sin7)-2P3(sin7)] (19) 

o 

C20 = -aia2/9^(4cos^-l), (20) 

y 

1 2 

C22 = -Q;iQ;2/3^sin^6' - -q;iQ;2/3^-P2(cos7). (21) 

o o 

The first set of coefficients, which do not contain a, capture most of the relevant physics. 
For the sake of completeness, the second set gives the a-dependent terms cn, C13, C20 and 
C22 which must be included for calculations in which ai and 0:2 cannot be neglected. 

From these expressions, it is evident that a good choice of coordinate system matters 
a great deal. Using the system discussed above (angles 6 and 7, separation r), leads to 
remarkably simple expressions for the coefficients c. Other systems we explored generally 
contained a lot of associated Legendre polynomials. It is easily seen in the above expression 
that taking the limit 6^0 reproduces Kaiser-Hamilton's plane-parallel result. 



3. Discussion 

We have presented a derivation of the redshift space correlation function between 
galaxies at two infinitesimal volume elements separated by an arbitrary angle, and have 
obtained very simple closed form analytic expressions for the correlations. These equations 
show that the effects of a finite angle are quite important. These are illustrated graphically 
in Figure 2 which indicates how the redshift distortions change as the angle between the 
two lines of sight is increased from to 120 degrees. 

Since most of the next generation redshift surveys are wide angle with the majority 
of pairs separated at angles greater than 10 degrees, the use of these expressions will 
dramatically increase the amount of information which can be extracted from the surveys 
resulting in more robust measurements and smaller shot noise contamination. Also, if the 
survey is not contiguous but consists of several slices like the Las Campanas Survey (Landy 
etal 1996), or the 2dF survey which will be formed from hundreds of pencilbeams, this 
approach can use all the data together to estimate the infall distortions. 



Another promising approach to study the large scale behavior of the power spectrum is 
with building up a Karhunen-Loeve basis (Vogeley and Szalay 1996). The first step in that 
method is to subdivide the survey into small cells (in redshift space) and then construct 
their correlation matrix. Most pairs of cells will have a large relative angle. Thus, with the 
expressions presented here, computing the correlation matrix is a trivial exercise. Once 
the appropriate basis has been created, one can use the Fisher matrix to select the set of 
eigenmodes most sensitive to /3, yielding an optimal estimation of the value of 13. 

In conclusion, we believe that our results allow a much more elegant treatment of the 
redshift distortions in a general geometry than before with the added benefit that all galaxy 
pairs may be utilized to construct a signal. 

AS would like to acknowledge the hospitality of the University of Tokyo, where this 
work has been completed. Useful discussions with Y. Suto are also acknowledged. 
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Fig. 1. — A graphic representation of our coordinate system. The angle between the two 
hnes of sight is 29. The main hne-of-sight is taken to be the axis bisecting the angle 29, 
since this is the symmetry axis for the two directions fi and Y2- The angle 7 is defined by 
the main line of sight and the vector r. 
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r sin 7 (h Mpc 



r sin 7 {h Mpc' 



Fig. 2. — Contour plots of the redshift-space correlation functions in polar coordinate system 
r, 7 for fixed finite angles 6. A standard CDM model with COBE normalization and Hq = 70 
km s^^ Mpc^^ is assumed. Plots (a), (b), (c) and (d) correspond to 6* = 0,30,45 and 60 
degrees. Solid and dashed lines indicate positive and negative values, respectively. Contour 
spacings are Alog^gl^l = 0.25. Solid thick lines represent the value 100, 10, 1 and 0.1. 
Dashed thick lines represent the value -1 and -0.1. Blank areas correspond to unphysical 
regions 7 > 6*. 



